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Abstract 

This paper is devoted to the presentation and the study of a new 
nonlinear subdivision scheme eliminating the Gibbs oscillations close 
to discontinuities. Its convergence, stability and order of approxima- 
tion are analyzed. It is proved that this scheme converges towards 
limit functions of Holder regularity index larger than 1.192. Numeri- 
cal estimates provide an Holder regularity index of 2.438. Up to our 
knowledge, this subdivision scheme is the first one that achieves simul- 
taneously the control of the Gibbs phenomenon and regularity index 
larger than 1 for its limit functions. 

Key Words. Nonlinear subdivision scheme, limit function, regularity, 
stability, Gibbs phenomenon. 
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1 Introduction 

Subdivision schemes are useful tools for generating smooth curves and sur- 
faces. For convergent schemes, starting from discrete sets of control points 
and using basic rules of low complexity, curves or surfaces can be obtained as 
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limits (called limit functions) of sequences of points generated by recursive 
applications of the subdivision scheme. 

A simple example of subdivision scheme is the family of interpolatory 
subdivision schemes, based on Lagrange's interpolation that have been de- 
rived and analyzed in [TT]. Other example is the family of spline subdivision 
schemes related to spline spaces |S]. 

The four-point interpolatory scheme [16], [15] is a convergent linear 
scheme of the first family, involving four-point stencils at each subdivision, 
for which the limit function is at least in the space 1 C . The Chaikin algo- 
rithm [7] is an example of spline subdivision scheme, with lower complexity 
as the previous example and converging towards C 2 ~ functions 2 . 

For applications, for instance to computer aided geometric design or 
image processing, complexity and convergence/regularity are not the only 
quality criterion. The order of approximation, that characterizes the preci- 
sion of the scheme, is an other important one. Moreover, oscillations that 
could occur in the limit function at the vicinity of strongly variating data 
(coming from the sampling of discontinuous functions) , called Gibbs oscil- 
lation, are really undesirable. 

In the last decade, various attempts to improve the properties of linear 
subdivision schemes, have lead to nonlinear subdivision schemes. For such 
schemes, the subdivision rules become data dependant; in addition to the 
previously defined criteria, one should add a stability property that ensures 
that the nonlinear scheme is linearly affected by perturbations of the data 
(for linear schemes, the stability is a direct consequence of the convergence). 

For nonlinear subdivision schemes, very few results concerning conver- 
gence or stability are available, see for instance [5], [9], [12], [21], [10] and 

A large family of nonlinear subdivision schemes to which belong the 
ENO, WENO or PPH schemes [9], [I] is made by the schemes constructed 
as a perturbation of the four-point linear interpolatory Lagrange scheme 
based on centered degree 3 polynomial interpolation. These schemes are in- 
terpolatory subdivision schemes (i.e. based on interpolation rules) and are 
constructed to avoid the Gibbs oscillations occurring classically for linear 
interpolatory schemes (see Figure [1]) . The schemes of this family are un- 
fortunately characterized by a low regularity of the limit functions of type 

^or < a < 1, / G C Q (R) iff / is bounded and 3C > such that Vx, y G R, \f{x) - 

f(y)\<C\x-y\ a 

For a > 1, / G C a (R) iff / (W) is bounded and G C {a - [a]) where [a] is the integer 

part of a 

2 c a ~ = {/ g c 8 ,v/3 < a} 



2 



C . Moreover, the ENO scheme is unstable. 

In p3], a new linear four-point subdivision scheme was presented. Its 
refinement rule is based on local cubic interpolation followed by a shift of 
1/4 or, in other words an evaluation at positions 1/4 and 3/4 rather then the 
standard evaluation at 1/2. This new scheme was shown to be convergent 
towards a C 2 curve. 

The aim of this paper is to analyze the scheme obtained using the same 
trick (shift of 1/4) for the PPH-type schemes [3] which are derived modi- 
fying the classical four point interpolatory subdivision scheme substituting 
the harmonic mean to the arithmetic mean. After the definition of the new 
scheme in section [2] we analyze successively its convergence (section [3]) , its 
stability (section 2]) and its order of approximation (in section [5]). Its be- 
havior in presence of strongly variating data (Gibbs oscillations) is analyzed 
in section [6l The last section is devoted to concluding remarks. 



2 A new nonlinear subdivision scheme 

As mentioned above, the starting point of our work is the construction of N. 
Dyn, M.S. Floater and K. Hormann in |14j . There, a new linear four-point 
subdivision scheme that generates C 2 curves is presented. Its refinement rule 
is based on the local cubic Lagrange interpolation, followed by evaluation at 
positions 1/4 and 3/4 of the refined interval. For all / € Z°°(Z), the scheme 
is then given by 

WJ)2n — 128 128 128 128 

(57)271+1 = ~i2gf n ~ 1 + i28^ n + l28^ n+1 ~~ 128^ n+2 ' ^ 

Following [4j where a nonlinear scheme is derived modifying the classical 
four-point interpolatory subdivision scheme substituting the harmonic mean 
to the arithmetic mean, we first obtain two new formulations of the scheme 



fQf , 49 f^ 14 f 1 f 7 (d 2 f n + d 2 f n+1 ) 

, cn 15 . 50 1 5 (d 2 f n + d 2 f n+1 ) 

WJ )2n+l — 77 In + 77/n+l _ rj/n+2 



64 1 "' 64 ~ 64 J " 64 
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2 

(Qf\ - f 15 f M 50 f 5 (d 2 fn + d 2 f n+1 ) 

PJj2n - 64 Jn-l + g4 Jn+ 64 Jn+l ^ 2 

fqf s _ 1. 49 f * 4 f 7 (d?f n + (Pf n+1 ) 

{^J)2n+l ~ 64 /"-l+ g 4 /n+ 64 /™+l g 4 g 

where (d 2 /) is defined by d 2 / n = / n+1 - 2/ n + / n _ x . 

The two formulations differs essentially in the distribution of the points 
f n contributing to the three first terms of 1 and 2. 

Using the same strategy as in [3], we define the new nonlinear subdivision 
scheme Sppha associated to (HJ by 

If \d 2 f n \ > IdVn+ll, 

(5'pPHA/)2n = g^/n + g-^/n+l + ~~ ^ PPH (d 2 fn, d 2 f n+ \) , 

(5'pPHA/)2n+l = g^/rt + - - ^ 4 PPH(d 2 /n,d 2 /n+l)- 



If \d 2 f n \ < \d 2 f n+1 \, 



1 . 15 „ . 50 „ 5 /j2 , j2 _ 



(^PPHA/)2n = —^fn-1 + g^/n + g^/n+l ~ ^ PPH ( d /«' d /n+0 

1 49 14 7 2 2 

(SpPHA/)2n+l = + g^/n + g^/n+l ~~ 64 PPH ^ ^ n '^ /i+0j 

where pph stands for the harmonic mean defined by 

(x, y) eM 2 ^ pph(x, y) := -^-(sgn(xy) + 1), (2) 

x + y 

with sgn(x) = 1 if x > and sgn(x) = — 1 if x < 0. 

The initial motivation for the substitution of the arithmetic mean by the 
harmonic mean is the elimination of oscillations near strong variating data 
thanks to the fact that 

|pPH(a:,y)| < 2min(|s|, \y\), (3) 

substitutes to 

x + y ^ n i i n 
— ^— < max{\x\, \y\). 

Before analyzing in details the properties of the new scheme 5"ppha 
we summarize the most important properties of the harmonic mean in the 
following proposition (see [5] for more details). 
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Proposition 1 For all (x,y) € M 2 , the harmonic mean PPn(x,y) satisfies 

1. ppn(x,y) = ppn(y,x). 

2. PPH(x,y) = ifxy<0. 

3. pph(— x, — y) = — pph(x, y). 



x-y 

x+y 



4. PPn(x,y) = Sl9n{x) + Si9n(y) min(\x\, \y\) [l + 

5. |pph(x, y)\ < max (\x\, \y\). 

6. \ppn(x,y)\ < 2min (\x\, \y\). 

7. For x,y > ; mm(x,y) < PPH(x,y) < ^5^. 

8. If x = 0(1), y = 0(1), \y - x\ = 0(h) and xy > then 

\^y-PPH(x,y)\=0(h 2 ). 

9. |ppH(xi,yi) - PPH(x 2 ,y 2 )| < 2max(|xi - x 2 \, \yi - y 2 |)- 
3 Convergence and Regularity 

We recall the following definition. 

Definition 1 A subdivision scheme S is said to be convergent if 

V/ G Z°°(Z), 3S°°f G C°(R) such that lim su Pn& \(S j f) n -S°° f(n2~ j ) \ = 0. 

(4) 

In order to derive the convergence, we rewrite the nonlinear subdivision 
scheme Sppha as a perturbation of a classical two-point linear subdivision 
scheme, 5c, introduced by G. Chaikin in [7j and defined by 



3 1 

(Scfhn = (5) 
1 3 

0Sc/)2n+l = 

The scheme Sc is known to be convergent with a regularity C 2 ~. 
Writing 



If \d 2 f n \ > \d l fn+l\, 



(S P PHA/)2n = ~Jn + ~rf n+1 + ^ fn+1 ~ TrPPH^/n,^/^), 
4 4 64 64 

13 1 5 

(£pPHA/)2n+l = ^/n + ~ ^d f n+1 - — PPH(d f n ,d f n +l), 

If \d 2 f n \ < \d 2 f n+1 \, 

(5 P PHA/)2n = -gfn + \fn+l ~ g^/n ~ ^ PPH ( rf2 /n, d 2 /n+l), 

13 1 7 

(SPPHA/Wl = ^/n + 4 + (54 ^V" ~ ^fn+l), 

we get that Sppha can be expressed as 

SpPHAf = S C f + F(d 2 f), 

with 

f m2« _ / md 2 f n+ i - ^PPH(d 2 f n ,d 2 f n+1 ) if \d 2 f n \ > \d 2 f n+ i\, 

( J)2n ~ I ~hd 2 f n+l - l 4 PPH(d 2 f n ,d 2 f n+1 ) if \d 2 f n \ < \d?f n+1 \, 

(6) 



and 



■2f\_ _ _ / -<k d2 fn+l ~ f|PPH(d 2 /n,d 2 /n+l) if \d 2 f n \ > \d 2 f n +l\, 

+ll- 
(7) 



^ /j2n+1 " 1 e^/n+l " l 4 PPK{d 2 f n ,d 2 f n ) if |d 2 /n| < |^/n + l|. 



To analyze the convergence of SpphAj we use a result proved in [3], [2] 
that reads: 

A sufficient condition for the convergence of a nonlinear subdivision 
scheme S NL : Z°°(Z) -> Z°°(Z) of the form: 

V/eP(Z), VnGZ ( (WWi=(S/) 2n+1 +F( ( 5/) 2n+1 , 

[ {jNLj)2n=Jn, 

where -F is a nonlinear operator defined on Z°°(Z), (5 is a linear and continuous 
operator on and S 1 is a linear and convergent subdivision scheme is: 

Theorem 1 If F, S and 5 given in (0j verify: 
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3M > such that \fd e Z°°(Z) ||-F(d)||oc < MHd^, (9) 
3c<l«icfrttorf ||55 , (/)+<5F(5/)|| 00 < cll^/lloo, (10) 

then the subdivision scheme Snl is uniformly convergent. Moreover, if S 
is C a convergent then, for all sequence f € i°°(Z), S^ L {f) G C' 3 wni/i 
/5 = min (a, —/052(c)). 

Using theorem [TJ we are going to prove the following result. 

Theorem 2 The nonlinear subdivision scheme SppHA * s convergent with a 
regularity at least C@~ with (3 > — log 2 {j^) > 1- 

Proof 

From the properties of the harmonic mean (Proposition [T]) , 

\dd 2 f n+ i - c 2 PPn(d 2 f n ,d 2 f n+1 )\ < max(ci,c 2 )||fi 2 /||oo. (11) 



For the perturbation F defined in © and ([7]), it is then easy to see that for 
alldeP°(Z), 

Halloo <^N|oo, (12) 

that is hypothesis ©. 

We now consider hypothesis (|10|) related, in this case, to the contraction 
of the second order differences (d 2 f). To simplify the notations we call 

f 1 = s PPHA (f). 

Different cases must be considered: 

Case 1: k=2n+l, study of (d 2 f v ) 2n+l = ft n+2 - 2fl n+1 + f\ n 

case lAi: \d 2 f n \ > \d 2 f n+1 \ and |c2 2 /„+ij > \d 2 f n +2\, 

case 1A 2 : \d 2 f n \ < \d 2 f n+1 \ and |d 2 /„+i| < \d 2 f n+2 \, 

case lBi: \d 2 f n \ > \d 2 f n+1 \ and |d 2 / n+ i| < \d 2 f n+2 \, 

case 1B 2 : \d 2 f n \ < |<i 2 /„ + i| and |d 2 / n+ i| > \d 2 f n+2 \. 



Case 2: k=2n, study of (d 2 f l ) 2n = ft n+1 - 2f\ n + f\ 



2n-l 



case 2Ai 
case 2A2 
case 2Bi 
case 2B2 



n I ■ 



\d 2 f n \ > \d 2 f n+1 \ and \d 2 f n ^\ > \d 2 f 
\d 2 f n \ < \d 2 fn+i\ and \d 2 f n -i\ < \d 2 f n \, 
\d 2 f n \ > \d 2 f n +i\ and \d 2 f n „i\ < \d 2 f n \, 
\d 2 f n \ < \d 2 fn+i\ and |d 2 /„_i| > \d 2 f n \. 
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The others cases follow by symmetry. 



• Cases 1A: We obtain for the case lAi 

21 1 2 1 1 2 ^ 2 

(d f )z n +l = - + ^/n + -^d / n+2 + —d f n+ i 

7 3 

- — PPH(d 2 f n+l ,d 2 f n+2 ) + — PPH(d 2 /n,d 2 /n+l) 

= ^d 2 f n+1 + —d 2 f n+2 + ^d 2 f n+1 
7 3 

-^PPH(d 2 /«+l,d 2 /n+2) + — PPH(d 2 /n,d 2 / n+ i) 

7 3 

- — PPH(d 2 / n+1 ,(i 2 / n+ 2) + — PPH((i 2 / n ,d 2 / n+ i). 

Using equation ([TTjl for d 2 f n+ i and PPH(d 2 / n+ i, d 2 f n+ 2) and Proposition [H 
we have 

K^Wll < 19+ 6 4 + 3 |l^l|oo<g||^/||oo. (13) 

Similarly for the case IA2, we have 

(d 2 / 1 ) 2 n+i = ^d 2 f n+1 - ^d 2 f n - ^PPR(d 2 f n+1 ,d 2 f n+2 ) 
g 

+—ppn(d 2 f n ,d 2 f n+1 ), 
o4 

and (1131) remains valid. 



• Cases IB: We obtain for the case lBi 

(d 2 f 1 )2n+l = ^/n+2 - + ^/n ~ ^d 2 /n+l + 64^ 2 ^ n+1 

5 3 

- — PPH((Z 2 /„ +1 ,d 2 / n+2 ) + — PPH(d 2 / n ,d 2 / n+ i) 

= ^d 2 f n+1 

5 3 
—^PPn(d 2 f n+l ,d 2 f n+2 ) + —ppH(d 2 f n ,d 2 f n+1 ). 
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Using equation (fTTj) for d 2 f n +i and PPH(<f 2 / n+ i, d 2 f n+ 2), and Proposi- 
tion [TJ we have 

K^/V+i! < -^piM 2 /lloc< ^llcP/Hoo. (14) 

Similarly for the case IB2, we have 

(d 2 f 1 )2n+l = ^^ 2 /n+l + -^d 2 f n+2 ~ g^^/fi 

-^PPH(d 2 /n+l,d 2 /n+ 2 ) + ^PPH(d 2 / n ,d 2 / n+ i), 



and (1141) remains valid. 



• Cases 2A: We obtain for the case 2Ai 

(d 2 f 1 )2n = 7/n+l - -Jn + 7/n-l _ ZT^/n+l - — d 2 f n 

4 4 4 o4 64 

9 5 
+ — PPH(d 2 / n ,d 2 / n+1 ) - — PPH(d 2 / n _i,d 2 / n ) 

= — fri — (i fn+] d fn 

9 5 

+ — PPH(d 2 f n ,d 2 f n+1 ) - — PPH((i 2 / n _i,d 2 / n ) 

= u d2fn ~ 64 d2fn+1 

+ -^PPH{d 2 f n ,d 2 f n+1 ) - ^PPH(d 2 / n _i,d 2 / n ). 

Using equation (fTTj) for c? 2 / n+ i and PPti(<i 2 / n , d 2 f n+ \), and for d 2 f n -\ and 
ppH(ii 2 / n _i,(i 2 / n ), we have 

\(d 2 f l )2n\ < ^^||rf 2 /||oo<^||rf 2 /||oo. (15) 

Similarly for the case 2A2, we have 
1 Q 1 

(d 2 / 1 )2n = f/fn + jr/fn-l 

3 7 

+ — PPH(d 2 f n ,d 2 f n+1 ) - —PPH(d 2 f n -l,d 2 f n ), 



and (1151) remains valid. 
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• Cases 2B: We obtain for the case 2Bi 

(d 2 f 1 )2n = -gfn+1 ~ -gfn + ~ g4^ 2 -^ n + 1 + ' fn-l 

9 7 

+ — PPH(d 2 f n ,d 2 f n+1 ) - — PPH((i 2 /n-l,d 2 /n) 

13 1 

= 4 d2fn ~ u d2fn+l + 64 d2/n_1 
9 7 

+ — PPH(d 2 f n ,d 2 f n+1 ) - — PPH((i 2 /n-l,d 2 /n). 

Using equation (TTTT) for d 2 f n+ i and PPn(d 2 f n ,d 2 f n+ i), and for d 2 / n and 
ppH(ii 2 / n _i,(i 2 / n ), we have 

l(d 2 / 1 )2n| < 16+ 6 4 + 1 |M 2 /lloc<g|M 2 /IU. (16) 
Similarly for the case 2B2, we have 

3 5 
+777 pp n( d2 fn,d 2 f n+1 ) - —PPK(d 2 f n - U d 2 f n ), 

04 o4 

and (|16j) remains valid. 

From equations (fl"3|) . (JHJ), (fl"5|) and (fl"6j) . we deduce that for all / £ 

Z°°(Z) 

||rf 2 5pPHA/||oo<^||d 2 /||oo. (17) 
lb 

Therefore Sppha verifies the hypothesis (jlOh of the Theorem [U In par- 
ticular, we obtain the convergence of Sppha- 

For the regularity, we use again Theorem [TJ According to the values a = 
2 and c = we obtain the regularity constant j3 = min (2,—log2 (jq)) ~ 
1.192. 

□ 

Numerical Regularity 

Following [20], the regularity of a limit function can be evaluated nu- 
merically. Using Si and S2 the subdivision schemes for the differences of 
order 1 and 2 associated to Sppha (that can be derived due to the specific 
definition of Sppha) , the following quantities are estimated for k = 1,2, 

-log 2 ( 2 fc[l ^ +1/)n+1 ~~ ( 5 fc +1 -^ n ll° l 



\\{S k f) n +i - {S 3 k f)n\\oo 
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They provide an estimate for (3\ and 02 such that the limit function belongs 
to C 1+/3l ~ and C 2+ ^ 2 ~. From table [TJ the numerical estimate of the regu- 
larity is C 2 ' 438- . We recall that the corresponding estimate for the linear 
scheme [H] is C 2 - 67 ~. 



3 


5 


6 


7 


8 


9 


10 


01 


0.9999 


0.9999 


1 


1 


1 


1 


02 


0.4395 


0.7738 


1.2615 


0.6541 


0.4387 


0.4388 



Table 1: Numerical estimates of the limit function regularity C l+ ^ and 

C 2 +P2 - for SppHA _ 



4 Stability 

For simplicity in notations we call for any initial sequence /°, and any j G N, 
f 3+ = S(fi). We recall the following definition. 

Definition 2 A convergent subdivision scheme is stable if 

3C <+oo such that Vf°,g° G \\S°°f - S°°g\\ L oa < C\\f° - g°\\i«>. 

(18) 

As for the convergence, to derive the stability of Sppha we use the 
following theorem of [2]. 

Theorem 3 If F,S and 5 given in (0) verify: 3M > 0, c < 1 such that 
Vf,9,di,d2, 

\\Fidx) - F(d 2 )\\ 00 < M\\d! - d 2 \U (19) 
\\S(S NL f - W)I|oq < c||<5(/ - g)\\oo, (20) 

then the nonlinear subdivision scheme Snl is stable. 

In order to check the hypotheses of Theorem [3] for Snl = <Sppha we 
first prove the following lemma. 

Lemma 1 Let be f, g G l°°(Z), if \d 2 f n \ > \d 2 f n+1 \ and \d 2 g n+1 \ > \d 2 g n \ 
then 

\d 2 f n+1 + d 2 g n - 2ppa( y d 2 g n ,d 2 g n+1 )\ < 3\\d 2 f - d 2 g\\ 00 . 
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Proof 

We consider different cases. 

If d 2 g n+ \d 2 g n < 0, we have 

d 2 f n+1 + d 2 g n - 2ppn(d 2 g n ,d 2 g n+1 ) = d 2 f n+1 + d 2 g n . 
* if d 2 f n+1 d 2 g n < 0, using that \d 2 f n \ > \d 2 f n+1 \ 

\d 2 f n +i + d 2 g n \ < \d 2 f n+ i - d 2 g n \ < \d 2 f n - d 2 g n \. 

2 fn+lC 

\d 2 g n \, we obtain 



* if d 2 f n j r id 2 g n > 0, we have d 2 f n +id 2 g n +i < and with \d 2 g n -\-i\ > 



\d 2 f n +i + d 2 g n \ < \d 2 f n+ i - d 2 g n+ i\. 



• If d 2 g n+ id 2 g n > 0, we recall (Proposition [[]) that if x, y > 0, min(x,?/) < 
PPH(x,y) < max(x,y). 

Without loss of generality, we suppose that d 2 g n > 0. We denote 



* if H > 0, 



H = d 2 f n+1 + d 2 g n - 2ppa(d 2 g n , d 2 g n+1 ). 

H < d 2 f n+ i + d 2 g n -2mm(\d 2 g n \,\d 2 g n+ i\) 

< d 2 f n+1 - d 2 g n 

< \d 2 f n -d 2 g n \. 



* if H < 0, 



H > d 2 f n+1 + d 2 g n - 2max (\d 2 g n \, \d 2 g n+1 \ 

> d 2 f n+ i + d 2 g n - 2d 2 g n+ i 

> d 2 f n+1 - d 2 g n+ i + d 2 g n - d 2 g n+1 . 



We have again to consider different cases according to the sign of d 2 f n . 

- d 2 f n > 0, we have d 2 f n - d 2 f n+ i > 

H > (d 2 f n+1 - d 2 g n+1 ) + (d 2 g n - d 2 g n+1 ) - (d 2 f n - d 2 f n+1 ) 
> 2(d 2 /n+i - d 2 g n+1 ) + (d 2 g n - d 2 f n ) > -3||d 2 / - d 2 g\\ 00 . 
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- d 2 f n < and d 2 f n+1 d 2 f n > 0, then d 2 f n+1 < 

H > (d 2 f n+1 - d 2 g n+ i) + (d 2 f n+1 - d 2 g n+1 ) 
> 2(d 2 f n+1 - d 2 g n+l ) > -2\\d 2 f - d 2 g\\ 00 . 

— d 2 f n < and d 2 f n +\d 2 f n < 0, we have from hypothesis that 

d 2 f n+1 + d 2 f n < 0, 

H > (d 2 f n+1 -d 2 g n+1 ) + (d° n -d 2 g n+1 )+d 2 f n+1 + d 2 f n 

> 2(d 2 f n+ i - d 2 g n+ {) + d 2 f n + d 2 g n 

> 2{d 2 f n+1 - d 2 g n+1 ) + d 2 f n - d 2 g n > -3\\d 2 f - d 2 g\\ 00 . 

□ 

We are now ready to prove the stability of 5ppha- 
Theorem 4 The scheme Sppha * s stable. 
Proof 

We check the hypotheses of Theorem [3l 

Firstly, we start with the hypothesis (|19p for F. 
Using the expressions of perturbation F, (|6|) and 0, and Proposition [H we 
obtain for all d±,d 2 G l°°(Z) that 

\\F(dx) - F(d 2 )\\oo < 1+ J' 2 |Mi-^lloo- 

64 

Secondly, we have to verify the contraction hypothesis (J20D - 

For a couple f,g€ P°(Z), we study {d 2 f 1 - d 2 g l )k for k = 2n + 1 (case 

1) or k = 2n (case 2). 

We consider different cases, according to the proof of Theorem [2] for / 

or g. 

For k=2n+l, we have 7 cases to study, see table [2j the others cases being 
deduced by symmetry. 



• Case 1A\ — 1A\: with equation (|13p we can obtain directly 

, , 2 ,i ,21 | . 19 + l + 7-2 + 3-2 [Ug , 2 
\d hn+l ~d g 2 n+i\ < \\df-dg 

< l\\d 2 f-d 2 g\\ O0 . 
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/ verifies 


g verifies 


notation 




/ verifies 


g verifies 


notation 




lAi 


lAi 


lAi - lAi 




IBi 


1A 2 


IBi - 1A 2 


Case 


lAi 


1A 2 


lAi - 1A 2 


Case 


IBi 


1B X 


IBi - XB± 


1A 


1A X 


IBi 


lAi - 15i 


IB 


IBi 


1B 2 


IB X - 1B 2 




lAi 


1B 2 


lAi - LB 2 




1B 2 


1A 2 


1B 2 - 1A 2 



Table 2: Cases to consider for k = 2n + 1 in the proof of the stability of 
SpPHA- 



Case 1A\ — 1A 2 : from equations (fT3"j) and (fT4"j) . we obtain 



d 2 f 2n +\ ~ d 2 g 2n +i 



19 



15 



64^ fn+1 ~ 64 d 9n+1 + M d fn+2 + 64^ ffr ' 



7 

64 



PPH(d 2 /n+l,C? 2 /n+2) 



5 3 9 

+— PPH(d 2 g n+1 ,d 2 g n+2 ) + — PPH(d 2 / n , d 2 f n+ i) + — pph^V, d 2 5n+ i) 

— {d 2 f n+ i - d 2 g n+l ) - —(ppn(d 2 f n+1 ,d 2 f n+2 ) - PPn(d 2 g n+1 ,d 2 g n+2 )) 

+^- A {d 2 f n+2 + d 2 g n+ i - 2ppH(d 2 g n+1 ,d 2 g n+2 )) 
64 
3 

+— (pPH(ci 2 / n ,d 2 / n+ i) - PPu(d 2 g n+1 ,d 2 g n+2 )) 
64 
3 

+T~Ad 2 9n + d 2 f n +i - 2ppn(d 2 g n , d 2 g n+1 )). 
64 



Applying Lemma Q] and Proposition [H we obtain 
\d 2 fL + i-d 2 g\ n+l \ < 



16 + 7-2 + 2 + l + 3-2 + 3-3 llj2j , j2 . 

\\d 2 f-d 2 g\ 



64 



- 64" J y| 



Similarly 

• Case lAi — 1B\: from equations (fl~3|) and (fl~4"J) . we obtain applying the 
Lemma Q] and the Proposition Q] 

\d 2 f l 2n+l -d 2 g\ n+1 \ < |l|d 2 /- Alloc- 



• Case lAi — 1B 2 : from equations (|13j) and (fT5j) . we obtain applying the 
Lemma [J and the Proposition [1] 

23 

\d 2 fL+i~ d 2 gl +1 \ < -\\d 2 f-d 2 g\\ OQ . 
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• Case lBi — IA2: from equations and dHJ), we obtain applying the 
Lemma Q] and the Proposition Q] 

\d 2 fl +1 -d 2 gl l+1 \ < |||rf 2 /-rf 2 5l|oo. 



Case lBi — 1B\: from equation ()14|) we conclude 

|d /an+i ~ d 92n+i\ < ^ \\df-dg\ 

< ^\\d 2 f-d 2 g\U 



• Case lBi — IB2: from equations (fl"4"j) and (fl"5j) . we obtain applying the 
Lemma [T]and the Proposition [1] 

72^1 j2„1 



• Case li?2 — I-A2 : from equations (|15p and (|14p . we obtain applying the 
Lemma Q] and the Proposition [1] 

l^ 2 /2n+l -d*92n+l\ < J£ 1 1^/ ~ d 2 g\ |oo- 



For k=2n, we have 7 other cases to study, see table El The others cases 
are deduced by symmetry. 

We notice that in equations (|14p and (|15p the same coefficients, but the 
cases are not completely equivalent. 





/ verifies 


5 verifies 


notation 




/ verifies 


g verifies 


notation 




2Ai 


2Ai 


2Ai - 2Ai 




2Bi 


2A 2 


2B 1 - 2A 2 


Case 


2Ai 


2A 2 


2A 1 - 2A 2 


Case 


2Bi 


2Bi 


2Bi - 2B X 


2A 


2Ai 


2Bi 


2Ai - 2Bi 


2B 


2Bi 


2B 2 


2Bi - 2B 2 




2Ai 


2B 2 


2A X - 1B 2 




2B 2 


2A 2 


2B 2 - 2A 2 



Table 3: Cases to consider for k = 2n in the proof of the stability of Sppha- 
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• Case 2A% — 2A\: with equation (|15p we can conclude 

ij2*i j2 l i 15 + 3 + 9 • 2 + 5 -2..- l2 

23 

< — IMV-^lloo- 

• Case 2^4i — 2^4 2 ■ from equations (|15|) and (|16j) . we obtain 

3 5 7 

-— PPH(dV, An+i) - ^rPPH(d 2 /n-i,d 2 /n) + PPH^Sn-l , d 2 g n )) 
64 64 64 

16 9 

= 27 (^V* _ d2 9n) + 7TT(PPH(d 2 /„,d 2 / n+ i) - PPH(d 2 5n , An-l)) 
64 64 
3 

- — (d 2 f n+1 + - 2ppH(d 2 s„, dV+i)) 
64 

5 

- — (pPH(d 2 /„_i,d 2 / n ) - PPH(d 2 g n -i,d 2 g n )) 
64 

-^-Ad 2 fn + d 2 g n ^i - 2ppn(d 2 g n ^ 1 ,d 2 g n )). 
64 

Applying the Lemma Q] and the Proposition [H we obtain 

t r2 j-1 ,2 1 I - 16 + 9 • 2 + 3 • 3 + 5 • 2 + 3 . . , 2 „ .9 I, 

\d 2 fL-d 2 gL\ < ^ IKY-c^lloo 

OQ 



Similarly, 

• Case 2^4i — 2i?i: from equations f) 151) and (|16p . we obtain applying the 
Lemma [1] and the Proposition [1] 

l^n-^Ll < Alloc- 



• Case 2.Ai — 2i?2: from equations (jl5[) and (fTTj) . we obtain applying the 
Lemma Q] and the Proposition Q] 

1 3 

\d 2 fL-d 2 g\ n \ < _(|d 2 /- Alloo- 
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• Case 2B\ — 2A 2 : from equations (fTB"|) and (fT6j) . we obtain applying the 
Lemma Q] and the Proposition [1] 

\d 2 fL-d 2 9l\ < f^f-d^U 



Case 2B\ — 2B\\ from equation (|16j) we conclude 

, -2 -1 ,21, . 16 + 3 + 1 + 9- 2 + 7- 2 ,,^^ 2 
\df2n-dg 2 n\ < \\df-dg\ 

< §||d 2 /-^l|oo- 



• Case 2B\ — 2B 2 : from equations (|16j) and (|T7l) . we obtain applying the 
Lemma Q] and the Proposition Q] 

l^n-^nl < ^\\d 2 !-d 2 g\U. 



• Case 2i?2 — 2A 2 : from equations (|17|) and (|16|) . we obtain applying the 
lemma Q] and the proposition [1] 

\d 2 fL-d 2 gl n \ < |||d 2 /-d 2 5l|oo- 



Finally, the hypotheses of theorem (|20p are verified and stability can be 
established. 
□ 

5 Order of approximation 

In this section we consider the reproduction of polynomials and the order of 
approximation of Sppha- 

We recall the following definitions. 

Definition 3 A subdivision scheme S is said to reproduce polynomials of 
degree k if for all polynomial P of degree k: 

if Vn E Z, f n = P(n), then 3P a polynomial of degree k such that (Sf) n = P(2~ l n). 
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Definition 4 A subdivision scheme S is said to have an order k of appro- 
ximation if for all function g G C k and all h > 0, 

if f = g(h.), then \Sf - g{2~ 1 h.)\ < Ch k . 

We then have the following property. 

Proposition 2 Sppha reproduces the polynomials of degree 2 with transla- 
tion of j. 

Proof 

We remark that for any P, polynomial of degree 2, and p = (P(n)) ne z, 
we have 

,J2 ,2 \ Pn+Pn+1 

Therefore, for the initial sequence p = (p n )n£Z, <Sppha(p) coincides with the 
application to p of the linear scheme |14| . In particular, the results of N. 
Dyn, M.S. Floater and K. Hormann [14] can be applied and the property of 
definition [3] is satisfied with P(.) = P(. - 1/4). 
□ 

Concerning the order of approximation the following proposition holds. 
Proposition 3 For all function g € C 4 ([0, 1]) and h > 0, if 

then 

if d 2 f n d 2 f n+ i > for all n£Z, then 

\\(Sp P HAf)n-g(2- 1 h(n - i))!^ = 0(/> 4 ), 

otherwise 

!IOSppHA/)n - g(2~ l h(n - i))!^ = 0(h 3 ). 

Proof 

According to Proposition [H we have that if d 2 f n d 2 f n+ \ > for all n G N 
then 

i „„/ ,2 t Jit \ d f n + d fn+1, /n /7,4\ 

|pph(c? f n ,d fn+i) 1 = 0{h ). 

Therefore, is S stands for the linear scheme defined in [14J, according to the 
definition of the .Sppha ; 



\\S PP EAf-Sf\\ 00 = 0(K 
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Since (see [2]) the scheme S is of order of approximation 4 we get the 
result when d f n d 2 fn+l > 0- I n other case, the reproduction of polynomials 
leads to 

||(5 P PHA/)n - g{2- l h{n - = 0(h 3 ). 

□ 

Remark 1 Following fWjj one can also establish, using the stability of .Sppha 
that \ \S^ pnA f-g\\ oo = 0(h 3 ). 



6 Elimination of the Gibbs phenomenon 

In this section we focus on the behavior of the scheme in presence of strongly 
variating data. The reference behavior deals with a step function as shown 
on Figure [TJ As it is visible on Figure [JJ left, high order linear schemes suffer 
from an oscillating behavior named as Gibbs phenomenon. 

According to D. Gottlieb and C.W. Shu [18], given a punctually discon- 
tinuous function / and its sampling f h defined by f% = f(nh), the Gibbs 
phenomenon deals with the convergence of S°°(f h ) towards /. It can be 
characterized by two features ([18] p. 244): 

1. Away from the discontinuity the convergence is rather slow and for 
any point x, 

\f(x)-S°°(f h )(x)\ = 0(h). 

2. There is an overshoot, close to the discontinuity, that does not diminish 
with reducing h; thus 

max \f{x) — S°°(f h )(x)\ does not tend to zero with h. 

We are now going to prove that the nonlinear schemes .Sppha does not 
suffer from the Gibbs phenomenon oscillations, as it can be guessed from 
Figure [TJ We have indeed the following 

Proposition 4 Given < £ < h, for any function f defined by: 

Vx < £, f{x) = f-(x) wiihf- e C°°(]-oc, £] , 
V* > e, f(x) = f+(x) withf + G C°°([e, +oo[ , 

and discontinuous in £, we have: 

• if \x\ > \K \f(x - |) - S™ PHA (f h )(x)\ = 0(h 3 ), 
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Figure 1: Comparison of limit functions for the same initial sequence (sam- 
pling of function (|2ip ). Left, linear scheme ([1]), right nonlinear scheme 
•Sppha 

• if \x\ < \K /-(0) + 0(h) > S? PHA (f h )(x) < f-(h) + 0(h). 
Proof 

Without loss of generality, we focus on [0, +oo[ and suppose that /_(£) > 

/+(£)• 

We first consider a single application of Sppha- Using Proposition [3] we 

get: 

. for n > 2 and n x G {2n, 2n + 1}, \S PPiiA (f h ) n i ~ f+i^Km - ±)| = 
0(h 3 ) 

• for n = 1 since / is discontinuous in £, d 2 f n = 0(1) and d 2 f n+ \ = 
0(h 2 ). Then, from Proposition [Q PPH(d 2 /n, d 2 f n+1 ) = 0(h 2 ). More- 
over, according to the definition of Sppha as a perturbation of the 
Chaikin scheme Sc we get that ISppHA^^ni — Sc(f h ) ni \ = 0(h 2 ). 
Since Sc is a second order scheme we get that \SppnA{f h )m — /+(2 _1 /i(rai 
I)| = 0(h 2 ) for nj G {2n, 2n + 1}. 

• for n = 0, d 2 f n d 2 f n+ \ < and therefore, according to the Proposition 
[Q ppn(d 2 f n , d 2 f n+ i) = 0. It is then easy to check, from the definition 
of Sppha that /-(0) < SppHA^n < SpPHAiWl < f+(h). However, 
writing Sppha as a perturbation of Sc we get that |<SppHA(/' l )ni — 
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Iterating, according to the stability of SppHA we get : 



• for x > \S\ 



PPHA 



'OO 



f\x)-f+(x- l/2)\ = 0(h 3 ). 



• for < x < |/t, the contraction of the second order differences (equa- 
tion (|17p ) and the fact that the Chaikin scheme Sc does not produce 
Gibbs oscillations allow to conclude. 

□ 

Before concluding this work, we come back to Figure Q] and to the com- 
parison between the limit functions obtained with 5"ppha and the limit func- 
tion obtained with linear subdivision schemes starting from the sampling f h 
of the discontinuous function: 



It appears from Figure [U that the nonlinear scheme Sppha exhibits a 
much better behavior close to the discontinuity than the linear scheme of 
comparable complexity. From Proposition [3] we know moreover that the 
scheme Sppha is of higher order than the Chaikin scheme. 

7 Conclusions 

In this paper, a new nonlinear subdivision scheme has been defined. It has 
many desirable properties. It is convergent with a regularity proved to be 
at least C 1192 ~ and numerically estimated at C 2,438_ . By construction, it 
is adapted to the presence of isolated discontinuities and the Gibbs phe- 
nomenon is eliminated. The scheme is also stable, that due the nonlinear 
nature is not a consequence of the convergence. Moreover, its order of con- 
vergence is 3. Recalling that it is constructed from a four-point centered 
stencil, all these properties make of this scheme a very good candidate for 
various applications. 
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